Fast mode of rotating atoms in one-dimensional lattice rings 
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We study the rotation of atoms in one-dimensional lattice rings. In particular, the "fast mode" , 
where the ground state atoms rotate faster than the stirring rotating the atoms, is studied both 
analytically and numerically. The conditions for the transition to the fast mode are found to be 
very different from that in continuum rings. We argue that these transition frequencies remain 
unchanged for bosonic condensates described in a mean field. We show that Fermionic interaction 
and filling factor have a significant effect on the transition to the fast mode, and Pauli principle may 
suppress it altogether. 



Recent progress in manipulating neutral atoms in- 
cludes paired states of fermions, the crossover between 
a Bose-Einstein condensate (BEC) molecules and a 
Bardeen-Cooper-SchriefTer (BCS) superfluid of loosely 
bound pairs P, S S 0, S B 0] ■ ^^.r, most of the exper- 
imental studies involving cold atoms were conducted in 
continuum. Optical lattices, however, allow many effects 
of interest: Superfluid to Mott insulator transition 3, 
Bloch oscillation of particles in lattices due to Bragg scat- 
tering 0, parametric atomic down conversion in BEC 
etc. 

By rotating atoms in a continuum, many phenom- 
ena have already been investigated: the appearance of 
vortices [H El in both the BEC and BCS 0, quan- 
tum Hall states for fermions with fast rotation fre- 
quency [IJand vortex lattices in the lowest Landau level 
for BECs|l^. [T^ . Experimentally, these rotations are re- 
alized by stirring the cold atoms. When the atoms are 
rotated in lattices, theoretical studies show that lattices 
lead to many new effects under rotation, such as struc- 
tural phase transitions of vortex matter [15,]. Also, near 
the superfluid-Mott insulator transition, the vortex core 
has a tendency toward the Mott insulating phase [l^ and 
second-order quantum phase transitions between states 
of different symmetries in a two dimensional (2D) lattice 
were observed at discrete rotation frequencies [itI . [Tsj . In 
particular, we have recently shown [1^ that it is possible 
for atoms to stay in a 2D lattice even if the rotation fre- 
quency is larger than the harmonic trapping frequency 
and density depletion in the trap can then be developed 
in such regime. 

However, in a 2D lattice, analytic results are not easy 
to obtain. The origin of some phenomena is thus hard 
to catch. The so-called fast mode for rotating atoms is 
such an exainple. This mode was shown to exist in 2D 
lattices [3; but we investigate this here in a ID 
lattice through both analysis and numerical simulations. 
Throughout this article, we define a current to be "posi- 
tive" if it is faster than the stirring force and "negative" 
if it is slower. Thus, the ground state (GS) is in the 
"fast mode" if its current is positive. While the ring lat- 
tice is much simpler, it catches much of the physics of 
the fast mode in higher dimensions, too. Specifically, we 



will study the origin of the fast mode, its dependence on 
quantum statistics and whether the fast mode is always 
present when the rotation frequency is large enough. Be- 
fore we move on, we remark that stirring a classical fluid 
can never achieve positive current in the rotating frame 
and thus the fast mode does not exist for a classical fluid. 

For completeness, we first review the rotation of a sin- 
gle quantum particle in a continuum ring with radius R. 
In the rotating frame, the Hamiltonian for such system 
is H ^ -2^^^ - ^Lzj where L^^ is the angular mo- 
mentum operator —ihd/d4>, m is the mass, (f) is the an- 
gular coordinate and R is the radius of the ring. H can 



be reduced to H 
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exp (incf) / \/2n , where n 



with eigenfunctions 4'n{<t') 

are integers. Except mVP'R^ /2 in H which obviously re- 
flects the centrifugal potential, the Hamiltonian is iden- 
tical to that of a q-charged particle put into a ring 
thread by a flux (j) = 2'KmflR^c/q The spectrum 
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shows that when m^lR^ /h is a half integer fc + 1/2, all en- 
ergy levels become two-fold degenerate, Ek+i+i — Ek-i; 
For all other values of fi, the ground states are not de- 
generate. Furthermore, when milR^/fi is a half inte- 
ger, the doubly degenerate ground states {n = 0, 1) can 
have either positive particle currents (as is the case for 
n = 1) or negative particle current (n — 0) in the ro- 
tating frame. For all other cases, the non-degenerate 
ground states {n — 0) have negative currents in the ro- 
tating frame. In what follows, we will show that lattices 
change many of these results. 

We now study the rotation of bosonic and fermionic 
atoms in a lattice ring using a Hubbard model. In the 
rotating frame, the single band Boson Hubbard Hamil- 
tonian is ^ITi, JJi, i21|] 
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H.c.) + Ub (1) 



where 6^ {hi) is the Boson creation (annihilation) op- 
erator at site i, and the modified single band Fermion 
Hubbard Hamiltonian is [13, [H, IS! 
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where c^^ {ci,a) is the creation (annihilation) operator for 
a fermion of spin a =t, i at site i. In the Eqs. |l]l2l), (i, j) 
indicates nearest neighbor pairs, t is the hopping term 
between neighboring sites, and H.c. means Hcrmitian 
conjugate. The interaction term for bosons (fcrmions) 

is c/b = UY.Mh -i)bth (Uf = c/Ei<T^*.T<iC^a) 

with U the interaction strength determined by the s-wave 
scattering length and the lattice potential. With equally 
spaced sites as exemplified by this paper, the geometric 
factor Kij is given hy K = f3 sin a/ 2, where a is the angle 
subtended by neighboring sites with respect to the axis 
of rotation, and /3 is a dimensionless constant of order 1 
characterizing the lattice geometry and depth (TtI. [TsI. [l9t . 
In Eqs. (I1I2|) . the lattice constant and h are set to be one, 
so energies are given in units of the hopping energy. Note 
that the Hamiltonian used in Eqs. (jll2p is an approxi- 
mation based on a perturbative treatment of fiL^: The 
Wannier basis states are the eigenstates of the Hamil- 
tonian without this term. While QLz may be included 
for an alternative (and possibly better) Wannier basis, 
Ref. (e.g., Fig. 3 therein) shows that the differences 
are negligible if is small. This, however, is the regime 
we are interested in here. 

As a building block for the following study, a sin- 
gle atom in the lattice with Na sites is studied first. 
Since wavefunctions have to be unique, the total wind- 
ing phase around the ring has to be a multiple of 
2tt, e.g., &{n) = n27r with n a nonnegative integer de- 
termining the symmetry of the state. The wavefunc- 
tions in an evenly distributed lattice site system are 
totally determined by the phase (j){n) = Q{n)/NA as 
IV^) = X]jexp(i(/)j)/V^ where |j) is the state 
when the particle is at site j. We can then use the 
wavefunctions to determine their energies, currents (see 
Fig. [T|) etc.. In particular, the energies for the low lying 
states are E{n) = —2 {tcos(j){n) + ilK sin (j){n)) , a linear 
function of rotation frequency for any fixed n. From 
E{n), we see that for Kil ^ t, the ground state has 
n — Tim = L-^a/4J, where [xj is the the largest integer 
number that is less than or equal to the specified number 
X. This means that any state with maximum phase differ- 
ences of 7r/2 between neighboring sites are ground state 
for large f2. Furthermore, the linear dependence of E{n) 
on Q should be contrasted to the quadratic dependence in 
continuum rings. Since the discrete rotational symmetry 
is not broken by the rotation (the fl term in Eqs. (dUl])), 
states with different symmetry experience level crossings 
as O changes. By comparing the energies of states with 
different n, the level crossings between these states can be 
determined analytically. Figure [1^ is an example of the 
energies and the level crossings for the low lying states 
of an atom in an 8-site ring lattice. These level crossings 
demonstrate how the system evolves between states of 
different rotational symmetries as V, increases [itI . [18| . 
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FIG. 1: Analytically calculated energies (a) currents (b) of 
the low lying states of an atom in the ring lattice. The inset 
in (a) shows a sketch of a rotating ring lattice with eight 
evenly distributed sites. "GJ" means the index of the low 
lying states. For example, GJ=1 is for the ground state and 
GJ=2 is for the first excited state, and so on. At Q > Qc, the 
GS current is a constant 2t, because the phase difi'erence n/2 
between neighboring-sites makes Q dependent term in Eq. Q 
have no contribution to the current. 

Corresponding to theses level crossings, the current 
and other observables change as well. As an important 
observable, the current is used in what follows to study 
the fast mode in the lattice ring. The currents are shown 
in Figure where the particle current between site i 
and j is calculated by fl?', Tsll 

Jij =i[n,j,_ffjj] ^ it{aia^ - H.c.) + KQ{aia^ + H.C.) (3) 

for each individual state. All the currents J plotted in 
this paper are the integrated currents along the ring. In 
a homogeneous ring, they are given by J = iVg Jij. It is 
not surprising to see the persistent currents of the excited 
states can be bigger than that of GS or have opposite di- 
rection. The excited states will be used later to construct 
the state of spin polarized fermions. 

The fast mode in the lattice ring happens when 
the winding number of GS is at its maximum nm, 
where neighboring sites have maximum phase difference 
n„i2T: / N A and thus maximum currents are achieved. By 
matching the energy of the state with winding number 
(n„i — 1) with the energy of the state with winding num- 
ber Urm we obtain Q.cK{l — cos27r/A^^) = t sm2'iT / N a to 
determine the rotational frequency above which GS is 
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in the fast mode. Setting Na = 4 reproduces the special 
result in Ref. [13) [3- From our general result, it can be 
seen that for larger Na, is higher to achieve the fast 
mode, which makes it difficult to be realized experimen- 
tally. In addition, a constant current of 2t is achieved 
when the phase difference of 7r/2 between neighboring 
sites makes the current independent of (see Eq. ^ and 
Fig.[T]). For comparison, particles in continuum rings can 
move faster than the stirring only when milB? /% is a half 
integer. The difference is because the lattice breaks the 
continuous rotational symmetry. 

A single atom in a few-site ring does not answer the 
question of whether the fast mode comes from finite num- 
ber of atoms or finite number of sites. Thus we study the 
same ring but with more than one atom, in which case the 
quantum statistics of the atoms comes into play. First 
we consider a BEC consisting of many zero-temperature 
bosons described by a mean field. When the ring is ho- 
mogeneous, the distribution of the atoms in the lattice 
does not depend on the rotation. This means that the 
mean field wavefunction is the same as that of the single- 
atom wavefunction except normalized to N atoms, and 
thus the transition to fast mode does not change. Since 
for a BEC many bosons are involved, the fast mode is 
therefore only due to finite number of sites in the lattice 
ring. 

Spin polarized fermions at zero temperature show more 
interesting dynamics. The GS current of N spin polarized 
fermions can be obtained by summing over the currents 
of the lowest N states of a single particle in the lattice, 
because these fermions, assuming no p-wave interaction, 
occupy the lowest N states according to Pauli's princi- 
ple. The resulting current is a linear function of rota- 
tion frequency and, at large can be either positive (eg. 
with 3 fermions) or negative (eg. with 2 fermions) as 
can be checked from Fig. ^p. Therefore, the fast mode 
may disappear even when there are only two spin polar- 
ized fermions. This disappearance of the fast mode is 
due to the Pauli principle, because Pauli principle pre- 
cludes occupying the same quantum state, including the 
ground state, by identical fermions. This will be further 
discussed in the following. 

When the fermions are not spin polarized, there can 
be s-wave interactions. The simplest example to include 
this interaction (see Eq. [2]) is to have one spin up fermion 
and one spin down fermion in the ring lattice. Fig. Oa. 
shows the current per fermion as a function of rotation 
frequency at various interactions for the two fermions. 
Similar to the single atom case, the current approaches 
constant 2t as the rotation frequency 17 increases. How- 
ever, this asymptotic process depends on the interaction 
and can be very slow. In general, the current is not nec- 
essarily a linear function of rotation frequency. It is im- 
portant to note that because of the increased dimension 
of the Hilbert space, there are more level crossings than 
that for single atom as shown in Fig. [21 Although there 
maybe accidental degeneracies, generally the rotation fre- 
quencies at which level crossings happen do not coincide 



for different U . At large rotation frequency (shown in 
Fig. [2b), the current is always positive, that is, the cur- 
rent changes continuously from noninteracting case to 
strongly interacting regime, to be compared with a sys- 
tem with more fermions discussed later in this paper. 
With even larger fi, the current will approach 2t for all 
U in this non-spin-polarized two-atom system. 
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FIG. 2: Numerically calculated current per fermion when one 
spin up fermion and one spin down fermion are rotated in 
an eight site-ring lattice. The current per atom always ap- 
proaches 2t when the rotation frequency is large enough. 

It is interesting to study combined effects of spin and 
interaction by putting four fermions with two up spins 
and two down spins in the lattice. Figure [3| shows the 
current as a function of the interaction at large rota- 
tion frequencies fi. With both strong attractive and re- 
pulsive interactions, positive currents and thus the fast 
mode persist. However, when the interaction strength 
is small J7 ~ 0, the currents surprisingly become nega- 
tive and are nearly linear in the rotation frequency. The 
case without interaction, ?7 = 0, is easy to understand, 
because the current is the sum over the currents of the 
lowest two states for spin up and spin down fermions 
shown in Fig. [TJd. From the discussion on the spin po- 
larized fermions, it is clear that the total current for the 
four noninteracting fermions is negative (thus not in fast 
mode) when is very large. Furthermore, Fig. [3] also 
shows that the absence of fast mode for non-spin polar- 
ized fermions extends to finite attractive and repulsive 
interactions. This means that the Pauli principle may 
suppress the fast mode in many-fermion system even at 
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nonzero interaction. 
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FIG. 3: Fast rotation of four fermions in 8 site-ring lattice. It 
shows numerically the effects of interaction and spin on the 
quantum phase transition between fast mode and non-fast 
mode. The non-fast mode around U ^ is caused by Fermi 
pressure (see text). 

Next, we try to explain the fast mode in such four- 
fermion system. At large attractive interaction ([/ <C 0), 
the fermions form local pairs (2^ . and thus behave like 
bosons, so the current is positive which is not surpris- 
ing given that fast rotating bosons always display a fast 
mode. Since the fermions get paired roughly indepen- 
dent of the rotation frequency, the transition from the 
fast mode to the non-fast mode happens at roughly the 
same U, independent of the rotating frequency as shown 
in Fig. [3] While with large repulsive interaction {U ^ 0), 
no such pairing exists and whether the transition happens 
at U depends on the rotation frequency. The larger Q is. 



the smaller repulsive U is required to make the transition 
happen (see Fig. [3]). Furthermore, since the fermions are 
not paired for strong repulsive interaction, the existence 
of fast modes at large rotation frequency means that the 
fast mode does not necessarily mean pairing or conden- 
sation of fermions. 

Although our study is meant to be an in-principle theo- 
retical discussion, a few comments about experiments are 
in order. The realization of a ring lattice with a tunable 
boundary phase twist has been proposed recently [23j . 
Since many bosons can be in the same fast mode, time 
of flight after turning off the potential of even a single 
lattice ring should allow detection of the momentum dis- 
tribution, and thus the fast mode, of the atoms. Detect- 
ing the fast mode of fermions may need additional twists, 
such as using multiple lattice rings, to enhance the sig- 
nals. Using multiple lattices to enhance the signal has 
been demonstrated recently 0|. 

To conclude, we have studied the rotation of bosonic 
and fermionic atoms in one dimensional lattice rings. We 
found that minimizing the ground state energy may give 
the fast mode and the transition to the fast mode in lat- 
tice rings is very different from that in continuum rings. 
Fermionic interaction and the filling factor are shown to 
have significant effects on the transitions to the fast mode 
and Fermi pressure may suppress the fast mode. Fi- 
nally, the fast mode is due to the finite number of lattice 
sites and is not associated to bosonic statistics, pairing 
of fermions, or superfluidity. 
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